Anderson localization in carbon nanotubes: defect density and temperature effects 
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The role of irradiation induced defects and temperature in the conducting properties of single- 
walled (10,10) carbon nanotubes has been analyzed by means of a first-principles approach. We find 
that di-vacancies modify strongly the energy dependence of the differential conductance, reducing 
also the number of contributing channels from two (ideal) to one. A small number of di-vacancies 
(5-9) brings up strong Anderson localization efi'ects and a seemly universal curve for the resistance 
as a function of the number of defects. It is also shown that low temperatures, around 15 — 65K, 
are enough to smooth out the fluctuations of the conductance without destroying the exponential 
dependence of the resistivity as a function of the tube length. 

PACS numbers: 73.63.Fg,73.23.-b,72.10.Fk 



The ubiquity of defects in materials science limits 

' somehow the performance of the material; this is even 
more dramatic when the dimensionality of the system 

', is reduced moving towards nanostructures, in particular 
carbon nanotubes Q . Understanding and controlling the 
conductance of these systems is decisive for their possible 
application in molecular devices 0, @, 0, tS;i ^ i 3.^ 9] • In 

\ perfect single- walled carbon nanotubes (SWNTs), elec- 
trons propagate ballistically if the inelastic processes can 
be neglected, i.e. if the electronic phase coherence length, 
Lrj,, is larger than the nanotube length, L. High quality 
metallic SWNTs exhibit as large as one micron • If 
inelastic processes, like the electron-phonon interaction, 

\ are important, the system conducts within the diffusive 
regime. On the other hand, it has been shown 0, |^ 
that defects induce Anderson localization in the electron 
states of a SWNT ii Lo < L < L^, Lq being the localiza- 
tion length. This is the regime we are interested in. 

This Anderson localization phenomenon is due to 
the concerted action o f q uantum confinement and de- 
fect scattering effects [lO|, . The theoretical calcu- 
lations presented in Ref. |8| were concentrated on the 
zero-temperature case. Temperature (T) adds new de- 
phasing mechanisms for the coherent electron transport 
in SWNTs that should be considered on the same foot- 
ing as the other sources of de-phasing. In this context, 
, a fully microscopic study of the electronic transport in 
defected SWNT in the Anderson localization regime tak- 
ing into account the effect of finite T is lacking. This 
is the goal of the present Letter: to fully understand 
the properties of the localization regime as a function 
of the nanowire length, temperature, density of defects 
and strength of the defect scattering potential for the 
paradigmatic (10,10)-chirality. 

As we are interested in modelling the role of defects 



in the conductivity of SWNTs, we have chosen to ana- 
lyze in detail the mono-vacancies and di-vacancies that 
are the most common defects that appear upon atom 
irradiation 9J; the comparison of these two cases will 
allow us to analyze the electron localization and its ef- 
fect on the nanotube conductance as a function of the 
different scattering potentials associated with those de- 
fects. As our calculations of the zero-bias conductance 
drop associated with a single defect show a conductance 
of 1.98Go for the mono- vacancy and 1.36Go for the lat- 
eral di-vacancy ( 2Go = 4e^/ft- is the conductance of an 
ideal nanotube), di-vacancies will show a more dramatic 
effect in the nanotube conductance Accordingly, al- 
though we will present results for mono-vacancies and 
di-vacancies in order to draw more general conclusions, 
most of our discussion will be concentrated on the di- 
vacancy case. 

The conductance calculations are done within a first- 
principles Local Orbital Density Functional (LODF) 
method that maps the full hamiltonian into a local- 
orbital one, allowing the use of the standard machinery 
of non-equilibrium Greens function in order to extract 
the transmission probability, and from that the nanotube 
conductance jl^ . The simulation geometry consists of a 
device region with the defected nanotubes connected to 
two semi-infinite perfect tubes. The geometry of the nan- 
otube around each vacancy is calculated using this LODF 
method and the defected nanotube is defined by adding 
N defects separated by a mean distance d. Then, the 
length of the nanotube, L, is L = N x d. The defects 
are assumed to be distributed along the rim of the tube 
with a distance between consecutive defects (measured 
along the tube axis) presenting a uniform random distri- 
bution between and 2d. The conductance calculations 
are performed for a statistically significant number of de- 
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feet realizations. For the sake of simplicity, we present re- 
sults for a nanotube with either only lateral di-vacancies 
(this is the most favorable geometry for di-vacancies 
or with only mono- vacancies. In our calculations, we ob- 
tain an effective local orbital hamiltonian associated with 
the sp^-basis set of the FIREBALL-orbitals details 
about the technique can be found in Ref.Q. Due to the 
local nature of this approach, it is feasible to calculate 
the electronic properties of very long nanotubes (up to 
several microns long) with an arbitrary distribution of 
defects. 

First we analyze the energy dependence of the nan- 
otube differential conductance, g{E), for two different 
density of defects and zero T. In Fig.l we plot g{E) of 
a (10,10) carbon nanotube for four random configura- 
tions of defects with two different inter-defect distances 
{d = 45.4nm and 75.5nm) and two number of defects (15 
and 25) introduced in the tube, typical values observed in 
the experiments In all cases we observe strong fluc- 
tuations of the differential conductance as a function of 
energy: this is a clear indication of the localized nature 
of the electronic states. Notice also, comparing panels 
(a) and (b) of Fig. 1, that g{E) also depends strongly on 
d. For finite T calculations the crucial point to realize 
is that electrons are injected into the nanotube with en- 
ergies within a window of the order of ksT. Then, it is 
expected that the strong fluctuations in the conductance 
would be washed out for a finite T. This will happen if the 
thermal coherent length Lt {Lt — fivp/ksT, vp being 
the Fermi velocity) is smaller than Lq. One naively say 
also that Anderson localization could not appear in this 
limit. This would be the case only if is of the order of 
Lo- As perfect metallic tubes exhibit ballistic response 
over long distances (microns), the electrons injected in 
the nanotube with different energies do not suffer inelas- 
tic scattering events and they behave as independent elec- 
trons |l4l |. In this picture, the total conductance is just 
the sum of many contributions from electrons injected 
within the thermal energy window |l5l| . 

In Fig. 2a we show the dependence of the nanotube 
resistance with the tube length for different temperatures 
and a given random distribution of di-vacancies. These 
curves have been calculated as follows. First, we evalu- 
ate the resistance of the longest tube (corresponding to 
iV = 25 defects in this case) with a given random config- 
uration of di-vacancies. Next, we calculate the resistance 
of a nanotube with N — 1 defects just eliminating the last 
di- vacancy (the length of the nanotube is then reduced to 
L ~ {N — l)d). This procedure is repeated until we end 
with a nanotube containing only one di-vacancy. Then, 
the whole calculation is repeated for several temperatures 
but with the same random configuration of defects. At 
zero T, as expected, the fluctuations in the resistance 
are very strong. The effect of temperature is striking: 
at room T, there is a complete reduction of the random 
fluctuations seen at T = OK. However we still get the 
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FIG. 1: (color on line) Differential conductance for two de- 
fected single-wall carbon nanotubes. Top panel, corresponds 
to a distance between defects of d = 45.4nm and the bot- 
tom panel for d — 75.5nm. Two different number of defects 
(15, continuous line, and 25, dashed line) are analyzed in both 
panels. Notice that the plotted energy window is smaller than 
the room temperature energy window (25meV). 



exponential increase of resistance with tube length that 
points out that the Anderson localization regime survives 
even at room T. Figure 2a also demonstrates that only a 
very low T (less than 20K) is needed to wash out the fluc- 
tuations of the resistance as a function of the tube length. 
This result explains why fluctuations are not seen at room 
T as experimentally reported in In Fig. 2b we plot 
both the zero and room T resistances for several random 
configurations of di- vacancies with d = 75.5nm. Notice 
that at room T the resistance of different configurations 
is very insensitive to the particular random distribution 
of di-vacancies. This is not the case for T = OK where 
there are strong fluctuations with respect to the average 
value. 

Figure 3 shows our calculations for the mean value of 
the room T resistance (as a result of an average over 
15 random cases, as explained above) as a function of 
the carbon nanotube length for different d's. Our results 
indicate that this room T resistance seems to fit a uni- 
versal curve once it is plotted in terms of the number of 
defects (iV) instead of the total length (L). An important 
result of performing a finite T calculation is that the re- 
sistance curves of Fig. 3 fit better as Rq xexp(L/Lo) than 
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FIG. 2; (color on line) a) Temperature dependence of the cal- 
culated resistance for a defected nanotube with d — 75.5nm. 
b) For T=0 (dashed lines) and T= 300K (continuous curves) 
we plot the calculated resistances for different random defect 
configurations with d = 75.5nm. The mean value of the resis- 
tances is also plotted in black. 



as i?o/2 X eyi-p{L/La) {Rq being 1/Go), for > 3 - 5. 
This is a strong indication of the reduction of the con- 
ducting channels induced by both defects and finite T 
El ll7|; only one channel is contributing to the total 
conductance. Moreover, this also suggests that the nan- 
otube enters in the localization regime for a very low 
density of defects, i.e. with more than 5-9 di-vacancies, 
provided the resulting Lq is smaller than L^. We have 
looked into the origin of the reduction in the number of 
channels in our system by diagonalizing the transmission 
matrix E3 used in our calculations. The diagonaliza- 
tion supports the previous discussion, although the full 
reduction to one-channel only appears once 5-9 defects 
have been introduced in the nanotube. Notice that Figs. 
1 and 2 already accounted for this result: in Fig.l, for 
both 25 or 15 defects the differential conductance is never 
larger than Gq; while in Fig. 2, the resistance in this case 
is never smaller than Rq for a number of defects larger 
than 5. 

The inset of Fig. 3 shows the effective localization 
length, Lqi as a function of d for di-vacancies at room T 
indicating that Lq ~ 8.5d. In Ref.f3| we found Lq « 4. Id 



after averaging the resistance in the log-scale over many 
different random cases; note that this is appropriate be- 
cause of the normal distribution of log{R) in ID sys- 
tems "^^^ localization length reflects, however, 
that at room T, and for a particular distribution of di- 
vacancies, the conductance of the system is calculated 
as the sum of the channel contributions associated with 
the energy window Ep ± fc^T (which represents a kind 
of effective average on energy of the conductance). This 
new estimation of Lq a-t room T slightly modify the av- 
erage distance between di-vacancies that were estimated 
in Ref. ;83 to have been created in the SWNTs by Ar+, 
suggesting that those values have to be reduced in all 
the cases by around 50%. Notice that this density of di- 
vacancies was obtained neglecting the effects that mono- 
vacancies are producing in the nanotube conductance. 
This is supported by our calculations for mono- vacancies 
whose results can be summarized in the following way : 
(i) the ballistic regime is found to extend up to around 
200 mono-vacancies ; (ii) the localization regime appear- 
ing for more than 200 defects shows a resistance as a 
function of the tube length that can also be fitted to 
i?o X ea;p(L/Lo) with Lq « 600d, implying that in the 
localization regime the conductance is again controlled 
by only one channel. Therefore, mono- vacancies show a 
much smaller efl'ect on the nanotube conductance than 
di-vacancies do. It is expected that for large vacancy 
clusters (larger than di-vacancies), Lq would be also pro- 
portional to d, but with a proportionality constant much 
smaller than 8.5. 
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FIG. 3: (color on line) Calculated room temperature resis- 
tance for different average distance between di-vacancies (d) 
as a function of number of defects. Inset: localization length, 
Z/o, extracted from the fitting R — B^y. exp(L/Lo) as a func- 
tion of d. 

Finally for the di-vacancy case we have analyzed the 
effect of temperature in the quenching of the conduc- 
tance fluctuations. This has been done for a tube with 
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25 di-vacancies as follows. As we are interested in the 
fluctuations that have already disappeared at room T, 
we take the calculated R{L) at room T as reference data. 
The scaled resistance plot is shown in the inset of Fig. 4 
for d = 75.5nm for a particular random configuration 
(the same as analyzed in Fig. 2a), after substraction of 
its mean value in L. We perform this calculation for all 
random configurations of defects and calculate for each 
of them the r.m.s deviation (cr) with respect to the R{L) 
data at room T. The final average over all configurations 
of fj is plotted in Fig. 4. The general trend is as follows: 
the higher T the faster the fluctuation damping. Also, 
the fluctuations decrease more rapidly for a low density 
of defects (large d). Figure 2a shows that for d = 75.5nm, 
the fluctuations are small for T higher than Tc ~ 15K, 
whereas for the other two d's (not shown here) this oc- 
curs for Tc « 32K and Tc « 65K, for d = 37.6nm and 
16.3nm, respectively. Thus, the temperature at which 
the fluctuations are quenched scales as 1/d. This can be 
understood by noting that at Tc the thermal length Lt 
should be of the order of the localization length, Lq (that 
is proportional to d) . Physically, the smaller the density 
of defects the smaller the averaged scattering potential, 
therefore a lower T would be needed to wash out the 
fluctuations. We have also carried out similar numerical 
calculations for (5,5) metallic SWNTs (not shown here). 
The underlying physics in this case seems to be the same 
than the one discussed for (10,10) SWNTs, allowing us 
to safely conclude that our theoretical findings would re- 
main valid for metallic SWNTs of different chirality. 



a higher number of defects the system shows localiza- 
tion, reducing the number of effective channels from two 
(ballistic) to one. (iii) At zero T, the conductance of the 
nanotube as a function of its length shows strong fluctua- 
tions. The net effect of a finite T is to wash out the strong 
fluctuations presented at T = OK; our calculations show 
that those fluctuations could be observed for di-vacancies 
(with d smaller than around 70nm and iV > 5 — 9) if the 
temperature is below 15 — 65K. We stress here that, in 
spite of the disappearance of the fluctuations, the expo- 
nential behavior of R{L) is still preserved at room T. This 
puts in evidence our assumption of a very low inelastic 
scattering in SWNTs, i.e. the phase coherence length of 
electrons much larger than the localization length. 

Our results are important in order to understand, in 
a further study, the high bias conductance where optical 
phonons do play a key role in limiting the conductance. 
Still there is a lack of knowledge about the role played 
by both temperature and defect density in this regime. 
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FIG. 4: (color on line) Calculated r.m.s deviation (a) as 
a function of temperature for three different density of di- 
vacancies. The inset shows the calculated resistance referred 
to the room temperature results for a particular random con- 
figuration with d — 75.5nm (see text for discussion). 

In conclusion, we have shown that in a (10,10) car- 
bon nanotube: (i) the transition between the ballistic 
and the localization regimes appears for a small number 
of di-vacancies in the nanotube (around 3 — 5). (ii) For 
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